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^^ ', Abstract 

1-^ I Recently, Benzoni-Gavage, Danchin, Descombes, and Jamet have 

Cd ' given a sufficient condition for linear and nonlinear stability of solitary 

wave solutions of Korteweg's model for phase-transitional isentropic 
gas dynamics in terms of convexity of a certain "moment of instability" 
I _^ ^ with respect to wave speed, which is equivalent to variational stability 

<^ • with respect to the associated Hamiltonian energy under a partial sub- 

1/-^ I set of the constraints of motion; they conjecture that this condition is 

t^^ ' also necessary. Here, we compute a sharp criterion for spectral stability 

f^ ' in terms of the second derivative of the Evans function at the origin, 

^. , o > 

^i ■ and show that it is equivalent to the variational condition obtained by 

vQ I Benzoni-Gavage et al, answering their conjecture in the positive. 

o' 
^. 

c^ . 1 Introduction 

Motivated by recent work of Benzoni-Gavage et al [BDD.Tj . we investigate 
in this paper stability of "soliton-type" (i.e., homoclinic) traveling wave 
/S ' solutions 

^ (1.1) U{x,t)=U{x-st), lim U{z) = Uoo, 

z— utoo 

U = {v,u), tJ = {v,u) of the Korteweg model 

vt — Ux = 

(1-2) * , ; 
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for isentropic phase-transitional gas dynamics, written in Lagrangian co- 
ordinates, with V denoting specific volume, u particle velocity, p pressure, 
and K > a coefficient of capillarity, taken for simplicity to be constant. 
The extension of our results to general k{v) as considered in |BDDJj is 
straightforward; see Remarks 11.51 and 11.71 As pointed out in |BSSj . system 
H1.2|) is also formally equivalent to the "good" Boussinesq equation modeling 
shallow-water flow. 

Equations ()1.1|) . accounting for the effects of compressibility and cap- 
illarity, but neglecting viscosity, are of dispersive, type, in contrast to the 
dissipative type of the usual compressible Navier-Stokes equations. Indeed, 
as discussed in greater generality in |BDD.Tj , (|l.lj) has the Hamiltonian struc- 
ture 

(1.3) Ut = JbU, 

where 

'0 1 



(1.4) J = 9,. J := 9,. , ^ Q 

is a constant-coefficient skew-symmetric first-order differential operator and 
bH, is a second-order differential operator corresponding to the variational 
derivative of the Hamiltonian functional 
(1-5) 

n= f H, H = (l/2)(n - uoo? - r {p{z) - p{voo))dz + {l/2)K{v){v^f, 

J Jvao 

of the (relative) total energy H of the system, with n{v) = k. constant. 
Formally, 

{d/dt)n{u) = {5n,Ut) = {5n,j6n) = o, 

so that the Hamiltonian is one conserved quantity of motion. A second 
(formally) conserved quantity, arising as a consequence of group invariance 
under translation (see |GSSj ) is the relative generalized momentum 

(1.6) Q{U) = {1/2){J{U - Uoo), {U - C/oo)) = j{u - u^){v - v^){x)dx 

(formally, J = {dxJ)^^dx, as prescribed in |(tSSJ ). Two additional conserved 
quantities are the relative masses 

Vi{U) = [{v - Vo^){x)dx, ViiU) = {u - Uoo){x)dx. 



The existence of these further quantities is associated with the fact that 
operator J is not onto, a circumstance that turns out to be significant. 

With this framework, it was shown by Benzoni-Gavage et al |BDDJj that 
stabihty of soUtons may be investigated by variational methods, fohowing 
the formahsm of jGSSl IBSj . Specifically, one may compute that solitary 
wave solutions are critical points of the Hamiltonian H, under constraint Q, 
satisfying the Euler-Lagrange equation (itself Hamiltonian) 

(1.7) {bn^sbQ){U) = ^, 

where speed s plays the role of a Lagrange multiplier. Such solutions occur 
in a one-parameter family IJ^ . Formally, strict variational stability of U with 
respect to constraint Q is thus sufficient for time- evolutionary orbital stabil- 
ity of the family {U^}, since, then, (i) A minimum U^^^i should therefore 
persist under small changes in Q, and (ii) within each level surface of Q, 
the conserved quantity W(C/)-7Y(0'^(2)) should control \\U -V''^^)f in the 
underlying Hilbert norm. 

In |GSSj . Grillakis et al made rigorous this intuitive argument, showing 
that variational stability indeed implies nonlinear orbital stability in a quite 
general framework (which applies here). Moreover, under the assumption 
that £* have at most one negative eigenvalue, they gave a necessary and 
sufficient condition for strict variational stability in terms of strict convexity 
with respect to s of the "moment of instability" 

(1.8) m{s) := {n + sQ){U') 
(in the terminology of |BSj ). i.e., 

(1.9) {(fm/ds^){s) = {d/ds)Q{W) = {{U - Uoo), J{dW/ds)) > 0. 

Finally, for J' onto, they showed by construction of a suitable Lyapunov 
function that strict failure of (|l.llj) implies time-evolutionary instability, 
completely deciding the issue of stability. 

Remark 1.1. Relation (|1.9|) shows that convexity of the moment of insta- 
bility is equivalent to monotone increase with respect to s of the generalized 
momentum Q along {U^}} Local monotonicity of Q along {tj^} is necessary 
for the picture of stability described above (in particular, for strict varia- 
tional stability), since the map from s to Q{U'^) must be locally invertible if 

^In the notation of |GSS| . monotone decrease with respect to lo = ~s. 



all sufficiently small perturbations (corresponding to small variations in Q) 
are to lie near some U^. Through the key relation 

(1.10) jC'idUyds) = -6Q 

(obtained by differentiating (|1.7|) with respect to s; see Proposition 4, |BDDJj ). 
we find that (|1.9I) is equivalent to 

(1.11) (fm/ds^ = -{dUyds,C'dUyds), 
where the self-adjoint operator £*, defined by 

(1.12) C' = {5^n + s6^Q)\us, 
is the constrained Hessian about U'^. 

In |BDDJj . a simple formula is given for m{s) and evaluated numerically 
to show that regions of both convexity and nonconvexity of m{-) may arise, 
depending on physical parameters. The first case corresponds with orbital 
stability, as discussed above. However, since J is not onto, the second case 
is inconclusive by the theory of |GSSj . The authors conjecture nonetheless 
that convexity is necessary as well as sufficient for stability, so that the 
second case in fact corresponds to instability. 

Here, we investigate time-evolutionary stability directly, using alterna- 
tive, Evans function methods introduced in [XTTjI IFWl IHZl IZll lB??Z] to 
obtain a sharp criterion for spectral time-evolutionary stability in terms of 
the sign of the second derivative at the origin of the Evans function D(X), 
an analytic function whose zeroes correspond in location and multiplicity 
with eigenvalues of L. Specifically, we obtain the following results deciding 
in the positive the conjecture of |BDD.Tj . 

Theorem 1.2. D"{0) = c{d?m/ds'^){s) for a nonzero constant c. 

Corollary 1.3. Traveling waves (|1.1() of H1.2() are time-evolutionarily (lin- 
early and nonlinearly) stable if d^m/ds"^ > and only if d'^m/ds'^ > 0, i.e., 
if they are strictly variationally stable with respect to constraint Q, and only 
if they are nonstrictly variationally stable. 

Theorem 11.21 generalizes a number of similar results obtained in |PWj 
for various related scalar models, and indeed holds in far greater generality 
for systems of abstract form H1.3|) . To see why, and to better understand in 



general the relation between the moment condition and variational and time- 
evolutionary stability, notice that, in this generality, 1)1.6(1 and (|1.9() become, 
formally (see |GSSj ^. 

Q = {l/2){J-^d,{U - C/oo), {U - Uoo) 

and 

icfm/ds^){s) = {j~^d,iu - f/oo), (duyds)) = {j-^d,u, (duyds)). 

On the other hand, differentiation of the traveling- wave equation with re- 
spect to X and s, respectively, reveals that /i = dxU is a right zero- 
eigenfunction of the linearized operator L and /2 = —dgU^ a generalized 
zero-eigenfunction of height two; for further discussion, see Section |^ Not- 
ing that left and right zero eigenfunctions /i and /i of L = J C are related, 
formally, by /i = J^^^fi, we find the general relation 

(1.13) {d^m/ds^){s) = -{hj2). 

But, vanishing of the inner product (|l.lc{|) of the genuine left eigenfunc- 
tion /i against the generalized right eigenfunction /i of height two precisely 
detects existence of a generalized eigenfunction of height three (by Jordan 
form), i.e., algebraic multiplicity of order two or more of the eigenvalue 
A = 0. Thus, {(Pm/ds^){s) must be a nonzero multiple of D"{0), since 
vanishing of D"{0) detects the same phenomenon. 

Formula (|1.13jl reveals a direct (formal) link between the moment con- 
dition and time-evolutionary stability that, moreover, does not require in- 
vertibility of J^ on the whole space, but only on the range of dx, through 
the distinguished roles of (dU^/dx) and —{dU^/ds) as zero eigenfunction 
and generalized eigenfunction of L. Recall that these same functions played 
critical roles in the argument of jCirSSj linking the moment condition and 
variational stability: (dU^/dx) as zero eigenfuction of C and —{dU^/ds) 
through relation (|1.10j) . 

This formal argument can be made rigorous using the extended spectral 
theory of |ZHj . Section 6, valid essentially wherever an Evans function can be 
analytically defined,^ concerning Jordan structure at eigenvalues embedded 
in essential spectrum, along with additional assumptions assuring that there 
exist no additional (extended) genuine zero eigenfunctions other than /o; 
see Remark 11.51 We shall instead follow the more concrete approach of 
direct Evans function calculations using the specific structure at hand, which 

^See also |MZIIZ2| for extensions to operators L of the degenerate type considered here. 



provide at the same time sign information. However, we note that these are 
quite similar to those on which the abstract development of |ZH| is based 
(see in particular the proof of Proposition 6.3, |ZHj ). 



Remark 1.4. Along similar lines, the general relation 

(1.14) 9fD(0)=c(/i,/,) 

has been established by Kapitula JK] for general (not necessarily Hamilto- 
nian) systems for which A = is an isolated eigenvalue of geometric multi- 
plicity one of the linearized operator L under consideration. In the present 
(Hamiltonian) case, A = is embedded in the essential spectrum of L, and 
so this result does not apply. Indeed, Benzoni-Gavage, Serre, and Zumbrun 
|BSZj have shown in the context of viscous conservation laws that (|1.14j) 
does not hold in general for embedded eigenvalue A = 0, but rather must be 
corrected by the addition of appropriate boundary terms at plus and minus 
spatial infinity; see also related results in KSl, KSSl lKKJ for perturbed NLS 
equations. Thus, the argument above reflects partly the special features of 
the Hamiltonian case. 

Remark 1.5. When ^7 is a differential operator of the general form dxJ 
considered here, with J an invertible (not necessarily differential) operator 
onto L^ (the most general case falling under the sufficient but not necessary 
theory of |GSSj l. we may take J~^dx = J~^ in the discussion above, and 

A = J-^fi = J-^dx{U - Uoo) = J~\U - f/oo). 

Noting that /i decays exponentially as x ^ ±00, we find that it is indeed a 
left genuine zero-eigenfunction of L, making rigorous sense of relation (|1.1,'-{|) . 
To complete the rest of the formal argument sketched above, note that 
bounded solutions of Lf = correspond to nearby homoclinic connections 
with different endstates, hence have equal limits as x — > ±00, and this 
eliminates them as possible extended eigenfunctions (which might in general 
be merely bounded |ZHj ) unless they in fact vanish at ±00. Thus, the 
standard assumption that [7 be a transverse connection, ensuring that the 
LP' kernel of L is one-dimensional, is sufficient to ensure that the extended 
kernel of L is also one-dimensional, i.e., there is a unique (up to constant 
factor) extended genuine eigenfunction /i with dual eigendirection /i. 

Moreover, since /i, /i, and /2 are all exponentially decaying, there 
exists an extended generalized eigenfunction /s, bounded but not necessarily 
decaying at infinity, if and only if (/i, 72) = 0. For, L/3 = /2 implies 

(/l,/2> = (/l,^/3) = (^7l,/3)=0. 



On the other hand, if the extended Jordan chain is order two, then (see 
|ZHj . Prop. 5.3 (iii)) the extended spectral projection 

^ = /l(/2,-)+/2(/l,-) 

preserves exponentiahy decaying elements of the extended eigenspace; in 
particular Vf2 = /2, whence (/i,/2) / 0. Finally, recalling QZHj . Theorem 
6.3 (ii)) that the Evans function vanishes with multiplicity equal to the 
dimension of the corresponding extended eigenspace, we obtain the result 

(1.15) D"{0) = c{<fm/ds^){s), c / 

for this general class: in particular, for the general isentropic Korteweg equa- 
tions discussed in |BDD,Tj consisting of (|1.3|) ~ (|1.5|) with arbitrary Hi{v) > 0. 

Remark 1.6. The same extended spectrum argument used in Remark ll.5l 
applied to general (not necessarily Hamiltonian) PDE, yields the remarkable 
fact that the dual version 

(1.16) dlD{0) = c{fk,h), 

c 7^ 0, of (|1.14j) remains valid in complete generality, without additional 
boundary terms, for embedded eigenvalues with (extended) geometric multi- 
plicity one for which the associated eigenfunction /i is exponentially decay- 
ing, a condition that is in the traveling wave context essentially always sat- 
isfied.'^ For, then Vfi = /i, where V = Ylj=i fjifK-j+i, •) and K > k is the 
order of vanishing of the Evans function, whence {fk, /i) = if and only if 
K > k+l, or, equivalently, d^D[0) = 0. Boundary terms arise in the forward 
formula through the integration by parts converting {fk,fi) = {fk,L^~^fk) 
to —l''~^{{L*)''~^fk,fk) = — l*'~^(/i,/fc), except in the case, as in Remark 
11.51 above, that fi, . . . , fk oi fi, . . . , fk decay exponentially. (Recall that ex- 
tended eigenfunctions fj, fj do not necessarily decay at infinity, except, by 
assumption, /i, but rather grow at most algebraically ZH , BSZ .) 

Remark 1.7. The relation (|1.15|) by itself does not give give the complete 
stability information of Corollary 11.31 but requires further information on 
the sign of c (with D suitably normalized: more precisely, the relation be- 
tween sgn(c) and sgnZ)(-(-oo)). In practice, this is often not restrictive, since 
we may calibrate the sign of c by continuation (homotopy) to a case where 

^Likewise, the obvious extensions to higher geometric multiplicity remain valid so long 
as all genuine extended eigenfunctions are exponentially decaying. 



stability is decided by the theory of |GSSj . i.e., either (^mjds^ > (vari- 
ational stability) or J onto. For example, in the case considered here, the 
collection of all soliton solutions of all isentropic Korteweg models comprises 
an open set in parameter space, so we may conclude by the (numerically ob- 
served) existence of stable waves. Alternatively, the regular perturbation 

J'^ := dxJ + OK, K := { n ) ' ^ sufficiently small, preserves sign infor- 

mation while converting J to an operator J that is onto for all 6^0. 

A still more general approach, not limited to the Hamiltonian setting, 
is to work directly from (jLlbf) . computing the sign by a direct calcula- 
tion like that of Section |31 In practice, one may sometimes determine the 
normalization without doing the full calculation, in which case it is easier 
to determine fk, /i numerically; see, for example, the analysis of stability 
of undercompressive traveling waves of thin- film models in jBMSZj (Prop. 
2.11 and footnote 6 on duality). We point out that the righthand side of 
H1.16() may by itself be considered as a generalized Melnikov integral, like 
(f'm/ds^ giving geometric information about the dynamics of the traveling- 
wave ODE; see |MZj Section 4.2 (especially eqns. (4.9)-(4.13)) for a general 
duality principle linking dual eigenfunctions to solutions of the adjoint ODE. 

Remark 1.8. There appears in |BDD,Tj the statement, apparently con- 
tradicting Corollarv 11.31 that there exist stable solitons that are variation- 
ally unstable. However, the time-evolutionarily stable waves considered in 
|BDDJj are in the standard sense variationally stable; indeed, this is the 
property that is used to prove time-evolutionary stability. The instability 
that is referred to, rather, is of the constrained Hessian, p.l2() . with respect 
to unconstrained perturbations, a notion that is necessary but not sufficient 
for instability with respect to constraint Q. 

Note. After the completion of this paper, we have learned of results of 
Bridges and Derks [BDl , BD2 establishing the relation D"{fi) = c{d'^m/ds'^) 
for the good Boussinesq system, which is equivalent to the main example 
H1.2|) studied here. More generally, they derive formulae generalizing those 
of Pego and Weinstein |PWj for scalar equations for the first nonvanishing 
derivative of a "symplectic Evans function" at the origin, for systems that 
can be put in the multi-symplectic form MZt + KZ^ = 'VS{Z), where 
Z G M^", M, K are constant skew-symmetric 2n x 2n matrices, and S" is a 
smooth function on M^". These include in principle a rather large class of 
Hamiltonian PDE, overlapping with but apparently distinct from the class 
(described in Remark II. 5 1) to which our methods apply. 



The multi-symplectic approach has the advantage that it is equation- 
independent once a change of coordinates to multi-symplectic form has been 
found, yielding automatically a characterization of the sign of c in (|1.15|) 
in terms of the geometry of the phase space of the traveling-wave ODE (cf. 
Remark ll.TL par. two). On the other hand, our methods are somewhat more 
straightforward, being carried out in the original coordinates and motivated 
by the simple relation (|l.i:^j) . In addition (see E,emarks ll.(im.7|) . they yield 
useful partial information also in the general, non-Hamiltonian case. 

2 Preliminaries 

Substituting Hl.l|) into H1.2() . we obtain the traveling- wave equation 

-sv' - n' = 
(2-1) . , y 

— SU + p[v) = —KV , 

or, substituting the first equation into the second, and integrating from — oo, 
the Hamiltonian ODE (nonlinear oscillator) 

(2.2) V" = K-l {s\ + p{v) - s\oc - p{Voo)) . 

Alternatively, we may write 1)2. 1|) formally as —sU' = JTC{U), or 

(2.3) J{6n + sJ){U) = J6{n + sQ){U) = 0. 

By the Hamiltonian structure of (|2.2|) . it follows that homoclinic orbits 
persist under changes in speed s and endstate C/oo forming for fixed C/qo 
a one-parameter family U^, s > 0, as described in the introduction. The 
equations are invariant under shifts in velocity u, so that any value of u^o 
is possible. However, the requirement that C/qo be a saddle-point of (|2.H) 
enforces on Voo the condition 

(2.4) s^+p'{voo)<0. 

Making the standard change of coordinates x — > x — si to a rest frame 
for the traveling wave, we may investigate its stability as an equilibrium 
solution U{x,t) = U{x) of Ut - sU^ = J6n{U), or 

(2.5) Ut = J6in + sQ){U). 
Linearizing (|2.5j) about U, we obtain 

(2.6) Ut = LU := JC, 

where £* is defined as in ()1.12|) and s denotes the speed of the wave U under 
investigation. 



Definition 2.1. The wave U is spectrally stable if the spectrum (t{L) of 
the linearized operator L about the wave is contained in the nonpositive 
complex half-plane {A : 3f?A < 0}. 

Routine calculation (see, e.g., |BDD.I1 RlT^Kl IFWJ . or Section El below) 
shows that the essential spectrum of L consists of the entire imaginary axis, 
so that stability is at best of neutral, or bounded type, rather than asymp- 
totic stability. This may be seen, likewise, by the fact that the equations are 
time-space reversible. Spectral stability is therefore determined by the point 
spectrum of L: specifically, whether there lie eigenvalues off of the imagi- 
nary axis. The following general result of Pego and Weinstein, a quantitative 
(linear) version of the previously-remarked relation between variational and 
time-evolutionary stability, gives a way to bound the number of such eigen- 
values. 



Lemma 2.2 ([PW]). For a linear operator L factoring as L = J L with 
J skew- symmetric and C self-adjoint, the number of eigenvalues of L in the 
positive complex half-plane {A : SRA > 0} is less than or equal to the number 
of negative eigenvalues of C. 

Proof See jPWl. D 

Lemma 12.21 was used in |BDD.Tj to establish the following upper bound. 
This is not needed in the present context, in which we seek to establish 
instability, but gives useful additional information; see Remark I.S.3I (It is 
necessary for stability; see Remark 12.41 ) 

Corollary 2.3 ([BDDJ]). The number of unstable (i.e., positive real part) 
eigenvalues of L = J Cf is less than or equal to one. 

Proof. Eigenvalues of C^ may be shown to correspond to eigenvalues of a 
(second-order scalar) Sturm-Liouville operator Mv := v" — k~^{s'^ + a)v, 
a := p'{v{x)) in the variable v. Since v' by translation invariance is a 
zero eigenfunction of this operator (see further discussion below; or, just 
differentiate (|2.1j) ). the number of unstable eigenvalues by standard Sturm- 
Liouville theory is equal to the number of nodes (zeroes) of v' , which, by 
(|2.1|) may be seen to be one. See BDDJ for details. D 



A general approach to resolving the issue of variational vs. evolutionary 
stability when variational methods fail to decide the question, introduced 
in |PWj . is via the Evans function D{X), an analytic function taking real 
values to real values, whose zeroes correspond to eigenvalues of L; for ori- 
gins of the Evans function, see |ET1 |E2l |EH| EH IA(t.TJ . By translational 
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invariance, -D(O) necessarily vanishes. For, differentiating (|2.3|) with respect 
to re, we obtain LU' = 0, hence U' is a zero eigenfunction of L. Likewise, 
existence of a one-parameter family [/** of solutions with the same endstate 
implies that -D'(O) = 0. For, differentiating H2.3() with respect to s yields 
L{dU^ /ds)\s=s = —U', hence —{dU^/ds)\s=s is a generalized zero eigenfmic- 
tion of L and there is a nontrivial Jordan block for L at A = 0. 

Since Z?(+oo) by standard considerations (see |A(t.T| IPW| KtZJ ) has a 
constant, nonzero sign as parameters are varied, this means that, in the 
generic case that D"{0) ^ 0, the number of unstable (i.e., positive) real 
roots of D (eigenvalues of L) is odd or even depending on the sign of D"{0). 
Used in conjunction with Lemma [2.21 this observation can yield complete in- 
formation. Namely, when the total number is at most one, spectral stability 
is determined by the sign of the second derivative of the Evans function. 

The second derivative D"{0) was evaluated in |PWj for scalar KdV- 
type equations and shown in several related cases to be exactly the second 
derivative with respect to s of the moment of instability, thus establishing 
necessity of (Pm/ds"^ > along with sufficiency oi(Pm/ds^ > for linearized 
and nonlinear stability in those cases. Their method of computing D and its 
derivatives does not apply in the system case considered here. However, it 
was shown in [£iZJ and jBSZj. respectively, that quite similar formulae may 
be obtained for -D'(O) and higher derivatives {d/d\)^D{Qi) in the system case; 



see also pfflTl IKK21 ITCm ITJ^ IRTTM IRD2J . 

In the remainder of the paper, we calculate D"{Q) by the method of 
|BSZj and show by explicit computation that it is a nonzero multiple of 
{d?m/ds'^){s){s), thus establishing Theorem 11.21 and Corollary 11.31 

Remark 2.4. The relation {d^m/ds^){s) = cD"{0), c / 0, shows that 
the assumption that £* have at most one negative eigenvalue is necessary 
for the result of |GSSj that (d'^m/ds'^){s) > implies stability. For, by the 
discussion above, the sign of I?"(0), hence of {d^m/ds^), counts the parity of 
the number of unstable eigenvalues, which might be even in general despite 
instability of L. 

3 Evans function calculations 

3.1 Construction of the Evans function 

Writing out the eigenvalue equation (L — A)[/ = 0, L as in (|2.()j) . we obtain 

\v — sv' — u = 
(3.1) 

Xu — su + (av) = —Kv , 

11 



where a := p'{v{x)) as in H2.1|) . which may be written in phase variables 
W = (n, V, v' ,v"y'^ as a first order system of ODE W = A{x, X)W, or 



(3.2) 



/n\' 



\v") 



( ° 


A 


-1 


\ 





1 














1 





V-A 


s\ 


—a — s 


-a'vxj 



/u\ 



V'J 



Using the fact (easily verified for p.l|) ) that U{x) converges exponentially to 
its limit Uoo as x ^ ±oo, we may construct an Evans function for (|3.2|) by 
the general method described in |GZj : see iZ2 for a more recent exposition. 
Examining first the constant-coefficient limit 



(3.3) 



Xv — sv' — u' = 
Xu — su' + acaV = —Kv'' 



of 1)3. If) . Ooo < by 1)2. 4|) . we find, taking the Fourier transform dx — > i^, 
that the spectrum of limiting, constant coefficient operator L^o of L satisfies 
dispersion relation 



(3.4) 



A(0 = siC ± Vooo - A^C^ e G 



hence consists of the imaginary axis. By a standard result of Henry |Hej on 
asymptotically constant-coefficient operators, we thus find that the essential 
spectrum of L likewise consists of the imaginary axis, while spectra lying 
in the strictly unstable half-plane A := {A : 5RA > 0} consists entirely of 
isolated eigenvalues. 

The same calculation shows that the dimensions of the stable and unsta- 
ble subspaces of the limiting coefficient matrix Aoo{X) := lim^^-i-oo ^(-z, A) 
is constant in A. For, substituting ^ = i^ into the characteristic equation 



(3.5) 



(A - s^if 



Q- 



-ooM^ + «^^^ 







determining eigenvalues /i(A), we obtain the dispersion relation (|3.4|) : thus, 
eigenvalues /u(A) may cross the imaginary axis only along the dispersion 
curve, i.e., for A on the imaginary axis. Taking A — > +oo along the real 
axis in (|3.5|) . we find that A^ ~ —Kfi'^, so that fi lie along the fourth roots 
of —1. Thus, both stable and unstable subspaces of Aoci{X) have dimension 
two in A. An elementary matrix perturbation calculation at A = shows 
that these extend analytically to a neighborhood of the origin, A = 0, with 
limiting values at A = given by the direct sum of stable (resp. unstable) 
subspace of A(X) and the vector 



(3.6) 



W = ( V^^, 1, 0, 0)*-- (resp. (- V=^, 1,0, 0) 



tr 
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with U = {v,u) coordinate equal to the unstable (resp. stable) subspace 

U = (1, V^^f' (resp. (1, -V^^)*' 
of the convection matrix 

-s -1 

(Recall, ()2.4() . that deto < 0, so that a has one real positive and one real 
negative eigenvalue.) 

Applying the framework of I GZj . we find that, on A, the subspaces 
of solutions decaying at +oo and — oo of the variable-coefficient equations 
H3.1|) likewise have dimension two. Moreover, there exist choices of bases 
W^ , W2+ and Wg", W^ for these subspaces that are analytic in A on A and 
extend analytically to A U B(0,r) for r > sufficiently small. At A = 0, 
we may choose W^{0,-) = 1^4(0, •)" = dxW^ asymptotically decaying, 
W^ := {u'^,v'^,v^,v^xY^, and W2^(0, •), W3"(0,-) asymptotically constant, 
with 

lim W^{x,0) = (c,l,0,0)*^ 

(3.7) "^+~ 

lim P^n-(x,0) = (-c,l,0,0)*^ 

I— >— 00 

c := \/—ctoo > 0. The Evans function is defined as the Wronskian 

(3.8) DiX) := det {W+ W+ W^ W7) |,=o, 

zeroes of which detect nontrivial intersection between decaying manifolds of 
solutions at plus/minus spatial infinity, i.e., decaying solutions of eigenvalue 
equation (|3.1j) . 

3.2 Calculation of D"(0) 

We now compute D"(0) by a simplified version (taking advantage of special 
structure) of the approach of |BSZj . 

Proof of Theorem M.'A Following |BSZj (see also related calculations of |ZHj , 
Section 6), we choose a convenient basis for the calculation of derivatives at 
the origin, based on the Jordan chain of L at A = 0. Namely, observing that 
Z := d\U satisfies at A = the first-order generalized eigenvalue equation 

(3.9) {L-\)Z = U 

13 



if U{X, •) satisfies for all A the eigenvalue equation {L — \)U = 0, we may 
arrange that not only H^-|^(0) = 1^4" (0) = dxW^ as above, but also 

(3.10) dxW^iO) = dxW^i^d) = -{dW'/ds){s). 

For, Wi and W^ may be chosen as "fast modes", decaying uniformly ex- 
ponentially as X — > ±00 for |A| sufficiently small, whence d\W^{^) and 
d\W^{^) are uniformly exponentially decaying at their associated spatial 
infinities as well, and (their U components) satisfy the generalized eigen- 
value equation (|3.9j) . which uniquely specifies them up to exponentially de- 
caying solutions of the eigenvalue equation: in this case, multiples ciVFj^ 
and C2W^ , respectively, which may be removed by the analytic change of 
coordinates Wj — > \cjWj. 

With this normalization, we find immediately that 

D{0) = det {d^W^ W^ Wg" d^W^) = 

and 

D'{0) = det {dxWj^ W^ Wr^ W^) +■■■ 
+ det{W+ W^ W^ d^W^) 
= dQi{dxW' W^ W3- dsW') 
+ det {dsW^ W^ W3- d^W') = 0. 

By a similar computation, we find that 

D"{0) = det {dlW^ W^ Wr^ W^) 

(3.12) +det{W^ W^ W3- dlW^) 

= det(l^+ W+ Wi {dlW^ - dlW+)) , 

W^ = dxW'^, where Yj := d\Uj satisfy at A = the second-order generalized 
eigenvalue equation 

(3.13) {L-\)Y = Z. 

(Recall that Z := d\U satisfies the first-order generalized eigenvalue equa- 
tion dSSl).) 

Now, setting A = and integrating (|3.1() from x = ±00 to x = 0, we 
obtain 

sv + u = isv + u)i^co), 

(3.14) ,, ^ '^ ' 
—su^av^KV = (— su -|- af)(±cxD), 
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for each {u,v) = (n, f)- associated with Wj . In particular, the righ- 
hand sides of ()3.14|) vanish for {u,v) = {u,v)f. Likewise, setting {u,v) = 
{u, v)- := d\{u, v)- , j = 1, 4, using (|3.1c{|) . and recalhng that we chose {u, v) 
to vanish at spatial infinity, we find that 

j-X 

sv + u = — / {dv^ /ds) dx, 



(3.15) 



-sii + av + Kv 



(du'/ds) dx. 



±oo 



Performing the row operations corresponding the the lefthand side of 
(|3.14j) . i.e., adding k~^{—su + av) to v", then adding sv to u, and using 
(Tmi) . (pmi . and (1X131) . reduces (1X1^ to 



(3.16) D"{0) = det 



/O s + c s — c Cu \ 

;+ f+ v^ {dlv^-dlvf) 



(vtY (4)' K)' (dlv^-dlvtr 
\ -sc - c^ sc- C^ Cu 



where Cv ■= — J^^idv^/ds) dx and Cu ■= n ^ J^^{du'^ /ds) dx. Since 



s + c s — c 

2 2 

-SC — C SC— C 



is invertible, by (|2.4|) . there exist constants d2, ds such that 

s + c 



^3 ( 2 

, SC — C 



do 



-SC — c 






Performing the corresponding column operation to eliminate C„, Cu, we 

obtain, finally, 

(3.17) 

/O s+cs-c \ 



D"{0) = det 



(^1+)' (4)' (^3")' (5>4- 



V 



2 2 

-SC — c sc — C 



dlvf + d2v:^ - d^v^y 



J 



Ci, 



where 

(3.18) 



C = det{ _l^_^^2 ,' _ '2 ) = Mp'i^oo) + s') < 0. 
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and 



(3.19) 7 - det ( ,^^^y (dlv^dlvt + d2vt - d,v7,)' ) ' 



Expand now 7 = 7_(0) — 7+(0), where 



(3.20) ,_(^):=det(j},, ^, 7+W - det ( Jj 



V{ Vl 



and 

(3.21) 7)4 := 9f-u^ + d3V^, Vl := dlvf + ^2^'^. 
Since ■()_ is bounded as x — > — cxd, while w^^ decays exponentially, 

(3.22) 7„(-cx))=0 
Similarly, 

(3.23) 7+(+oo) = 0. 
By (jXHl), (|XTT)|) . W+ = Ua; satisfies 

t;" + K~^(s^ + a)7; = 0, 



or 



v\' f 1\ fv\ fO 



^^•^^^ \v'J \-K-^{s^ + a) Oj\v'J~\0 

while v± satisfy 

v" + K-\s^ + a)v = K-^ I [-s{dv'/ds) + {du'/ds)] dx + C, 



or 



v\' / 1\ [ v\ /O 



^^•^^^ V^V V-«"'(s' + a) {))\v')~\F)' 

with C constant, 

F := K-^ I [-s{dv'/ds) + {du'/ds)] dx + C. 
J —00 
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Thus, using inhomogeneous Abel's formula and H3.22() - ()3.23|) . we may 
evaluate the difference 7 between Wronskians 7+ and 7^ at x = as the 
Melnikov-type integral 



(3.26) 



7 = 7-(0)-7+(0) 

F^ 

K^' I v^{x) f [-s{dv'/ds) + {du'/ds)]{y)dydx 
J —00 J —00 

K-^ / {v- Voo)[-s{dv'/ds) + {du'/ds)]{x) dx 



^detf'^r t\ix)dx 

-00 

/■+00 
-1 



00 



where in the second-to-last line we are using the fact that term C J_ Vx {x) dx 
integrates to zero and in the last line integration by parts. 

Substituting the relation —s{v — Uoo) = {u — Uoo) coming from the first 
equation in 1)2. 1() and recalling ()1.9|) . Remark II. 11 we obtain, finally, 

/+00 
[(n - Uoo){dvyds) + {v - Vco){duyds)]{x) dx 
-00 

^^•^'^^ =K-^ f °°(C7 - Uoo)J{dU'/ds){x) dx 



' —00 

-l/j2_ / j„2\ 



= K-'{d^m/ds'){s). 
Combining (|3.17j) and (|3.27|) . we obtain 
(3.28) D"{0) = {-C/K){d^m/ds^){s), 

—C/k > 0, completing the proof. D 

3.3 Proof of Corollary IT31 

By standard considerations |PW| I(tZ| IBSZ| IZ2j we have also the following. 

Lemma 3.1. As X ^ +00 along the real axis, sgnZ)(A) has limit +1. 

Proof. First, note that D{X) does not vanish for JRA sufficiently large, a stan- 
dard fact associated with well-posedness of the linearized time-evolutionary 
system; this may be established either by asymptotic ODE theory as in 
|(tZ| IZ2j or by elementary energy estimates as in [BSZ] . Thus, D{X) has a 
(nonzero) limiting sign as claimed. 
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To determine the value of this sign, we may consider a homotopy from 
system (|3.1() to the constant-coefficient equation 

Xv — u = 
(3-29) ^ ^ n> n 

Xu + KV = 0, 

capturing high-frequency behavior, for which bases 

Vj = {X/fij , 1 , ^j , /i^) , fij = 9j \/A 

of exponential solutions may be explicitly calculated for all A, where 9j are 
the fourth roots of —1. 

Choosing bases Vi, (V2 — Vi)/vA and V3, (V4 — Vaj/vA of decaying sub- 
spaces at ±00, we find that these extend continuously to A = 0, with the 
projection onto {u, v) coordinates of the limiting bases equal to the standard 
Euclidean basis (1,0), (0, 1); in particular, the projection of the limiting sub- 
spaces is nonsingular. Likewise, it is nonsingular for the bases at A = for 
the original system, and (by identical calculation) for all systems strictly be- 
tween the two. Thus, choosing continuously initializing bases at A = for the 
family of systems, we find that the signs of D{+oo) and of the determinants 
d± of the projections of decaying subspaces at ±00 onto {u, v) coordinates 
are all constant under homotopy, with, moreover sgnD(+oo)d-d+ explicitly 
calculable from system 1)3. 29() . We may therefore determine sgnZ?(+oo) for 
the original system by the straightforward calculation of d_ and fi+ , yielding 
the result. D 

Remark 3.2. An alternative, somewhat simpler approach, using the ho- 
moclinic structure of the wave, is to note that the limiting subspaces at plus 
and minus infinity of decaying solutions of the eigenvalue equation are com- 
plementary subspaces of the limiting coefficient matrix A, hence transverse. 
Thus, their determinant (real-valued, by construction) is of constant sign in- 
dependent of A, which may be related to sgnZ)(+oo) by explicit calculation. 
Calculating the value at A = of this determinant then gives the result. 
This is similar in spirit to the original argument of PW . However, this 
approach does not extend to the heteroclinic case, and so we have chosen to 
give the more general argument above. For related arguments, see |Z21IZ3j . 

We may now easily obtain the main result. 

Proof of Corollary ] 1. 'A It has been shown by variational considerations in 
|BDDJj that [d'^m/ds^){s) > implies linearized and nonlinear stability. 
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using Corollary 12.31 and the argument of |GSSj . Thus, we need only estab- 
lish that {(Pm/ds'^){s) < implies instability, for which it is sufficient to 
show that there is an L^ eigenfunction of L with positive real eigenvalue A. 
Defining the stability index 

(3.30) r := sgnL>"(0)i:>(+oo), 

we have, provided T ^ 0, that the number of positive real eigenvalues of 
L (zeroes of D) is odd or even, according as T is negative or positive, by 
standard degree theory on the line (recall, D is real-valued, by construction, 
for A real). Thus, F < implies existence of at least one positive real eigen- 
value, for which (see Section ^ the associated eigenfunction is necessarily 
exponentially decaying as x — > ±00. Noting that sgnF = sgn(d^m/ds^)(s) 
by (|3.28|) combined with Lemma l3. 11 we are done. D 

Remark 3.3. By Lemma 12.21 there is in fact precisely one unstable eigen- 
value, necessarily real, in case {d'^m/ds'^){s) < 0. 

4 Remarks on the viscous case 

Finally, we comment briefly on the viscous case, as modeled by 

vt — Ux = 
(4.1) * ; 

Ut + P[V)x = (-Uxx - KVxxx, 

e > 0. Stability of traveling waves of this model has been studied in detail 
in |ZlllOZj for a parameter range of (k, e) for which 1)4. 1|) may be converted 
by a change of dependent variables to a strictly parabolic problem. 

Including viscosity of arbitrarily small strength e > 0, one finds by the 
energy argument of (OZj that all homoclinic connections of the e = equa- 
tions (2) break, save for zero-speed solitary waves identical to those of ()1.2() . 
Moreover, by a calculation similar to but much simpler than that of the pre- 
vious section, we find that all of these waves are unstable, with -D(O) = 0, 
Z)'(O) = ce J{vx)'^dx, c 7^ 0. That is, even infinitesimal viscosity will either 
break or destabilize any solitary wave, whether stable or unstable for the 
e = model, leaving unclear the physical implications of stability or insta- 
bility with respect to (|1.2|) . However, depending whether the e = version 
is stable or unstable, the unstable root will be near or far from the origin, 
corresponding perhaps to some type of metastable phenomenon. 
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